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BACKGROUND OF THE INVENTION 



1 . Field of the Invention 

[0001] The present invention is related to the field of oil well logging. More 
5 specifically, it relates to improved techniques for data acquisition in induced gamma ray 
logging methods. 

2. Description of the Related Art 

[0002] Well logging systems have been utilized in hydrocarbon exploration for many 
years. Such systems provide data for use by geologists and petroleum engineers in 

10 making many determinations pertinent to hydrocarbon exploration. In particular, these 
systems provide data for subsurface structural mapping, defining the lithology of 
subsurface formations, identifying hydrocarbon-productive zones and interpreting 
reservoir characteristics and contents. Many types of well logging systems exist which 
measure different formation parameters such as conductivity, travel time of acoustic 

1 5 waves within the formation and the like. 

[0003] One class of well logging systems seeks to measure the incidence of nuclear 
particles on the well logging tool from the formation for analysis purposes well known in 
the art. These systems take various forms, including those measuring natural gamma rays 
from the formation. Still other systems measure gamma rays from the formation induced 
20 by bursts of neutrons entering the formation from a neutron source carried by the tool and 
pulsed at a preselected interval. 
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[0004] Additionally, it is frequently useful to obtain data regarding the time spectral 
distributions of the occurrence of the gamma rays. Such data can yield extremely 
valuable information about the formation, such as identification of lithologies which are 
potentially hydrocarbon-producing. Moreover, this desired spectral data may not only be 
5 limited to that of natural gamma rays, but also may be desired for the gamma ray spectra 
caused by bombardment of the formation with the aforementioned pulsed neutron 
sources. 

[0005] Well logging systems for measuring neutron absorption in a formation use a 
pulsed neutron source providing bursts of very fast, high-energy neutrons. Pulsing the 

10 neutron source permits the measurement of the macroscopic thermal neutron absorption 
capture cross-section E of a formation. The capture cross-section of a reservoir rock is 
indicative of the porosity, formation water salinity and the quantity and type of 
hydrocarbons contained in the pore spaces. In wireline well logging operations, as the 
neutron source is pulsed and the measurements made, the instrument is continuously 

15 pulled up through the borehole. This movement makes it possible to evaluate formation 
characteristics over a range of depths. 

[0006] The measurement of neutron population decay rate is made cyclically. The 
neutron source is pulsed for 40-60 microseconds to create a neutron population. 
Neutrons leaving the pulsed source interact with the surrounding environment and are 
20 slowed. In downhole environment, collisions between the neutrons and the surrounding 
fluid and formation atoms slow the neutrons. Such collisions may impart sufficient 
energy to these atoms to leave them in an excited state, from which gamma rays are 
emitted after a short time as the atoms return to a stable state. Such emitted gamma rays 
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are called inelastic gamma rays. As the neutrons are slowed to the thermal state, they 
may be captured by atoms in the surrounding matter. Atoms capturing such neutrons are 
also elevated to an excited state, and emit gamma rays after a short time as the atoms 
return to a stable state. Gamma rays emitted due to this neutron capture reaction are 
5 called capture gamma rays. 

[0007] Several prior art references are aimed at improving measurements in capture 
cross-section logging. For example, a method and apparatus employing a source and two 
detectors are discussed in U.S. Patent Nos. 4,645,926 and 4,656,354, both issued to 
Randall. A subsurface instrument includes a long-spaced (LS) and short-spaced (SS) 

10 detector for detecting natural or induced gamma ray emissions from subsurface 

formations. The detectors produce electrical pulses, each pulse corresponding in time 
with the incidence of a corresponding gamma ray on the detector and having an analog 
voltage amplitude correlative of the gamma ray. A method is discussed in Randall '354 
for determining presence of a gas by comparing first and second parameters obtained at 

15 the detectors. The first parameter is indicative of a count of detected impingements of 
primarily inelastic gamma radiation upon a detector. The second parameter is indicative 
of a count of detected impingements of primarily capture gamma radiation upon a 
detector. Randall '926 discusses a method of determining a parameter of the borehole, 
where primarily inelastic gamma radiation is normalized upon impingement on a 

20 detection device. 

[0008] The number of capture gamma rays present at any time is in direct proportion 
to the number of thermal neutrons, i.e., the thermal neutron population. The decay rate of 
this neutron population is an exponential function and may be defined by specifying the 
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time required for the thermal neutron population to decrease to one-half ("the 'half-life 
time'")- Although the neutron lifetime may be measured, a more useful parameter is the 
capture cross-section. Capture cross-section and neutron lifetime are inversely related, 
the capture cross-section being a measure of the rate at which thermal neutrons are 
5 captured in the formation. Analysis of a formation using the capture cross-section or the 
half-life time is referred to as "neutron decay analysis." 

[0009] The well logging instrument measures the capture gamma rays to determine 
the thermal neutron decay time or "neutron lifetime." The neutron lifetime represents 
combined effects of the thermal neutron capture cross-section in each of several volumes 
10 of space surrounding the logging instrument ("regions"). These regions generally 
include the instrument itself, the fluid in the wellbore, the steel casing, the cement, the 
earth formations radially proximal to the wellbore wall (which may have been infiltrated 
by fluid from within the wellbore), and the earth formations radially more distal from the 
wellbore wall (which have minimal infiltration from the fluid in the wellbore). 

15 [0010] In U.S. Patent No. 4,668,863, issued to Gray, et. al., an apparatus is used to 
analyze and process parameters including the macroscopic thermal neutron absorption 
capture cross-section of the formation at borehole elevations corresponding to the 
locations from which spectra are derived. For acquiring temporal spectral data, a multi- 
channel scale section is provided which includes a channel number generator which 

20 produces a numerical sequence of memory address codes corresponding to a sequence of 
adjacent time windows. Each code uniquely defines a start time. The windows 
collectively comprise the time interval of the desired spectrum. Each time a gamma ray 
pulse is detected, the tool addresses a memory location corresponding to the detection 
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time and increments the count value resident therein. At the conclusion of the time 
spectrum interval of interest, the memory locations may be interrogated by the CPU and 
the resultant spectral data analyzed, transmitted to the surface or presented visually as a 
gamma ray emission count versus time plot. Correlation is made of detection signals in 
5 response to impingement of gamma radiation upon first and second detectors. 

[0011 J A method employed by Randall et. al. in "Time-Derived Sigma for Pulsed 
Neutron Capture Logging" SPE 9614, p. 345, 1981, determines the timing of the measure 
interval using on a two-pass E measurement. That is, two E computations are 
continuously derived. A first-pass E value is measured over the fixed-time interval from 

10 400 /isec to 1000 /xsec after each neutron burst. Simultaneously, a second-pass E 
measurement is made over a single gate of 600 /zsec gate width which starts at a 
selectable time between 200 /xsec and 400 fisec after the burst. The positioning of the 
second-pass gate is dependent upon the first-pass E value. Although positioning an 
acquisition window can remove the effect of noisy data points from subsequent 

15 calculation, significant noise can still be introduced. 

[0012] There is a need for a technique of adjusting an acquisition window size to 
improve the quality of data used in gamma ray calculations. The present invention 
fulfills this need. 

SUMMARY OF THE INVENTION 

20 [0013] The present invention is a method and apparatus for determining a 

parameter of interest of an earth formation with a logging tool having a pulsed source 
used to irradiate the earth formation. A detector spaced apart from the source the results 
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of interaction of said irradiating with said earth formation. The received signals are 
analyzed using a processor that defines a starting time for analysis at a time at which non- 
formation effects are small and an ending time for analysis at a time at which the signals 
are statistically significant. 

5 [0014] In one embodiment of the invention, the source is a pulsed neutron source 

and a cross-section of the formation is determined. The starting time for analysis is 
based in part on signal levels towards the end of the analysis window corresponding to an 
earlier pulsing of the source. The ending time for analysis may be determined by 
comparing the signal level with a running average of the signals. 

10 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] Fig. 1 (prior art) shows a schematic of an induced gamma ray well logging 
instrument suitable for use with the present invention; 
Fig. 2 shows pulsed neutron decay data substantially free from noise; 
15 Fig. 3 shows noisy pulsed neutron decay data with an adjusted acquisition window; 
Fig. 4 shows noisy pulsed neutron decay data with an acquisition window of varying 
duration; and 

Fig. 5 shows a flowchart of a method of the present invention. 

DESCRIPTION OF PREFERRED EMBODIMENT 

20 [0016] Referring now to Fig 1, there is illustrated a nuclear well logging apparatus 
suitable for use with one embodiment of the present invention. This particular 
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configuration is an example of prior art that could be used and is not intended to be a 
limitation, and there are other arrangements that could be used. For example, many of 
the functions that are described herein as being done by electronic circuitry could also be 
performed by a suitable microprocessor. 

5 [0017] Well 10, also called a borehole, penetrates the earth's surface and may or may 
not be cased depending upon the particular well being investigated. Disposed within well 
10 is subsurface well logging instrument 12. The system diagramed in Fig. 1 is a 
microprocessor-based nuclear well logging system using multi-channel scale analysis for 
determining the timing distributions of the detected gamma rays. Well logging 

10 instrument 12 includes long-spaced (LS) detector 14, short-spaced (SS) detector 16 and 
pulsed neutron source 18. In one embodiment of the invention, LS and SS detectors 14 
and 16 are comprised of bismuth-germanate (BGO) crystals coupled to photomultiplier 
tubes. In another embodiment of the invention, the detectors are made of sodium iodide 
(Nal). Other suitable materials for detectors include Gadolinium Orthosilicate. To 

15 protect the detector systems from the high temperatures encountered in boreholes, the 
detector system may be mounted in a Dewar-type flask. Source 18 comprises a pulsed 
neutron source using a D-T reaction where deuterium ions are accelerated into a tritium 
target, thereby generating neutrons having an energy of approximately 14 MeV. Power is 
supplied to the neutron source 18 through power supply 15. Cable 20 suspends 

20 instrument 12 in well 10 and comprises conductors for electrically connecting instrument 
12 with the surface apparatus discussed below. Instead of an electrical connection, a 
fiber optic connection may be used for connecting the instrument 12 to the surface 
apparatus. 
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[0018] The outputs from LS and SS detectors 14 and 16 are coupled to detector board 
22 which amplifies these outputs and compares them to an adjustable discriminator level 
for passage to channel generator 26. Channel generator 26 is a component of multi- 
channel scale (MCS) section 24 which further includes spectrum accumulator 28 and 
5 downhole processor 30. MCS section 24 accumulates spectral data in spectrum 

accumulator 28 using a channel number generated by channel generator 26 and. After all 
of the channels have had their data accumulated, CPU 30 reads the spectrum, or 
collection of data from all of the channels, and sends the data to modem 32 which is 
coupled to cable 20 for transmission of the data over a communication link to the surface 
10 apparatus. Channel generator 26 generates synchronization signals which control the 
pulse frequency of source 18. The downhole processor 30 further communicates control 
commands which define certain operational parameters of instrument 12, including the 
discriminator levels of detector board 22, and the filament current and accelerator voltage 
supplied to source 18 by power supply 15. 

15 [0019] The surface apparatus includes may include a processor 33 coupled to cable 
20 for recovery of data from instrument 12 and for transmitting command signals to 
instrument 12. The surface apparatus also includes a depth controller 36 which provides 
signals to the processor 33 indicating the movement of instrument 12 within well 10. A 
display unit , plotter, and mass storage unit (not shown) may also be associated with the 

20 processor 34. The primary purpose of display unit and plotter is to provide visual 

indications of the generated logging data as well as systems operations data. Storage unit 
is provided for storing logging data generated by the system as well as for retrieval of 
stored data and system operation programs. The surface equipment may also include a 
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transceiver (not shown) for a satellite link enabling transmission of data to and receiving 
instructions from a remote location. 



[0020] In a well logging system such as is illustrated by Fig. 1, processor 33 initially 
transmits system operation programs and command signals to be implemented by 
5 downhole processor 30, such programs and signals relating to the particular well logging 
operation. Instrument 12 is then caused to traverse well 10 in a conventional manner, 
with source 18 being pulsed in response to synchronization signals from channel 
generator 26. Typically, source 18 is pulsed at a rate of 1000 bursts/second (1 KHz). 
This, in turn, causes a burst of high energy neutrons on the order of 14 MeV to be 

10 introduced into the surrounding formation to be investigated. In a manner previously 
described, this population of high energy neutrons introduced into the formation will 
cause the generation of gamma rays within the formation which at various times will 
impinge on LS and SS detectors 14 and 16. As each gamma ray thus impinges upon the 
crystal-photomultiplier tube arrangement of the detectors, a voltage pulse having an 

1 5 amplitude functionally related to the energy of the particular gamma ray is delivered to 
detector board 22. Dynamic temperature compensation is provided and a digital signal is 
passed to channel generator 26 of MCS section 24. 

[0021] Figure 2 shows typical pulsed neutron decay data for a formation sigma of 20 
20 c.u. and a borehole sigma of 100 c.u. Time is displayed in microseconds along the 

abscissa, and the count rate is displayed logarithmically in counts/second (cps) along the 
ordinate. Data are shown over a time interval of 0 - 1000 /is following a pulsing of the 
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source. From this data, a narrower processing window 203 is defined having a start time 
206 of approximately 370 /xs and an end time of approximately 970/zs. The data located 
within the processing window of Fig. 2 is substantially free from noise. Therefore, 
adjusting the acquisition window to improve accuracy is unnecessary in this case. These 
5 data can be reasonably processed with a full 600-/xsec wide processing window as 
indicated. As is the usual practice in the field of nuclear logging, the ordinate has a 
logarithmic scale. 

[0022] Figure 3 shows typical pulsed neutron decay data obtained from a wellbore 
comprising a formation having a E of 55 c.u. and a borehole having a E of 100 c.u. If the 

10 data of Fig. 3 are processed using a full 600-/xsec wide acquisition window beginning and 
ending at the times such as that indicated in Fig. 2, noisy data would be included, thereby 
reducing the accuracy of later calculations. Prior art methods that use a fixed processing 
window length would correct for the noise at the end of the acquisition by simply starting 
the processing window at an earlier time, such as 305 and end the processing window at 

15 306 and keeping a fixed processing window length of 600 ^s. Such prior art methods for 
positioning that acquisition window are discussed, for example, in Randall et al. "Time 
Derived Sigma for Pulsed Neutron Capture Logging", SPE 9614, , 1981. By acquiring 
data earlier in the cycle, some noisy data found late in the cycle can be disregarded. 
Thus, by having an earlier start time for the processing window, the noisy data acquired 

20 during the time period from 800 /isec to 1000 /xsec is not. Some noisy data is still 

obtained in the approximate time period of 700 /xsec to 800 /xsec. Moving the start of the 
processing window to still earlier times while avoiding late, noisy data, would lead to 
another problem. Data acquired at earlier times (before approximately 170 /xs ) is 
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typically affected by nuclei from the borehole as these nuclei have a faster decay rate. 
Therefore, moving a processing window to earlier times in order to avoid noisy data late 
in the cycle can lead to processing undesired data from the borehole that is not indicative 
of formation properties.. 

5 [0023] Figure 4 shows the same data of Figure 3 processed using a variable 

acquisition window width using the method of the present invention. As in Fig. 3, the 
processing window starts at an earlier time, -170 /xsec (405). The end of the processing 
window occurs at -700 /xsec (406). Thus, the duration of the processing window of Fig. 
4 is reduced (530 /xsec) over the window of Fig. 3. As a result, the noisy data in the 
10 region of 700 /xsec to 800 jxsec is removed from further calculations, thereby enabling 
improved accuracy, while, at the same time, avoiding undesirable signals from the 
borehole. Furthermore, the reduced processing window size of the present invention 
enables an operator to avoid undesired data points at both early and late times. 

[0024] The present invention combines previous methods with a statistical test to 
15 determine when the signal has dropped into the noise as described above and then 
narrows the processing window so that only meaningful data are used in the 
determination of the apparent formation E . For nuclear logging, the entire acquisition 
(and processing) window is described in terms of channels that correspond to short time 
windows within the acquisition window. A typical decay rate can be determined using 
20 the following equation: 
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iend 



i-istr 



0) 



iend 



where t is the decay rate, istr and iend are the start and end channels used for the data 
acquisition, tchfij is the time in the middle of the i-th channel i, cntfij is the number of 
5 counts recorded in the I th channel. In the present invention, the start channel (or the 
beginning of the processing window) for data acquired for the second and subsequent 
pulsing of the source is determined by using a sliding gate formula An earlier previous 
pulsing of the source (referred to as a sweep) gives rise to a late gate E value from the last 
channels of the sweep. In the current sweep, this late gate sigma value of the previous 
10 sweep is used to determine the time of the start of the processing window, using Eq. (2): 



In order to determine the time at which the end channel of this processing window 
occurs, an exponential running average is calculated and tracked. The exponential 
running average is an initialized sum average (sumavg = 0.0). As each consecutive 
15 channel is accessed, the average is modified in the appropriate manner. A typical 
modification scheme of the present invention is shown in Eq. (3): 



istr = 



2*4545.4545 



(2) 



late _ gate _ sigma 



sumavg = 0.8 * sumavg + 0.2 * cowitef/] 



(3) 



[0025] Figure 5 shows a flowchart of a method of the present invention. An 
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exponential average, sumend, of the last five channels of the previous sweep is computed 
in Box 501 and is used to determine a starting point for tracking the exponential running 
average of the current cycle. It is important to define the starting time for analysis so that 
the signals are not affected by non-formation contributions, such as from the borehole. 
5 An exponential average of the counts [i] is maintained beginning at the channel where the 
number of counts in the channel falls below twice sumend {counts fi]<2.0*sumend). This 
is shown in Box 503. The factor of two used in this step has been found to give the best 
results, though other factors could also be used. This average is referred to herein as 
sumy. Sumy represents the exponential running average from the present channel i 

10 (determined using counts fi]<2.0*sumend) to channel /+4. Sumy is updated as each 

consecutive channel is accumulated (Box 505). Once a channel is accessed in which the 
number of counts in the current channel registers below 0.8 of the current sumy 
(counts [i] <0.8*sumy) y as in Box 507, the spectrum is assumed to have fallen below a 
threshold at which the background noise masks any useful information. If the counts 

15 remains above this value the next channel is recorded (Box 506), and sumy is updated 
(Box 505). Otherwise the current channel becomes the end channel (Box 509). 

[0026] The processing window length is calculated from the algorithm of Eq. (4): 

iwidth = i - istr 

if (iwidth < 30), iwidth = (iwidth + 30)/2 (4) 
if (iwidth < 20), iwidth = 20 

The calculated length (iwidth) is then compared with the length from the previous scan 
20 (prevwidth). If iwidth > prevjwidth, iwidth = 0.75*prev_width+0.25*iwidth. Each 
calculation of iwidth is rounded up to the next integer. 
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[0027] Upon completion of data acquisition, the decay rate can be computed from Eq. 
(1). This decay rate t obtained is used to calculate the apparent formation E using Eq. (5): 

_ a x + a 2 i + a 3 t 2 + a/ 

L - " (5) 

where a/, a^, aj, and a 4 are functions of the channel width. The cross section is a 
5 parameter of interest that can be determined using the method of the present invention, 
invention. 

[0028] The invention has been described above with reference to an apparatus 
conveyed on a wireline. The method is equally useful when conveyed on coiled tubing 
for use in substantially horizontal boreholes, or on drilling tubulars (including drillstrings 
10 and coiled tubing) in a Measurement while Drilling (MWD) implementation. 

[0029] While the foregoing disclosure is directed to the preferred embodiments of the 
invention, various modifications will be apparent to those skilled in the art. It is intended 
that all such variations within the scope and spirit of the appended claims be embraced by 
the foregoing disclosure. 

15 
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